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Abstract 

Dirac's method for constraints is used for solving the problem of exclusion of double occupancy 
for Correlated Electrons. The constraints are enforced by the pair operator Q{x) = 'il:\^{x)ip^{x) 
which annihilates the ground state |^''^ >. Away from half fillings the operator Q{x) is replaced by a 
set of first class Non-Abelian constraints \x) restricted to negative energies. The propagator 
for a single hole away from half fillings is determined by modified measure which is a function of 
the time duration of the hole propagator. As a result: a) The imaginary part of the self energy - is 
linear in the frequency. At large hole concentrations a Fermi Liquid self energy is obtained, b) For 
the Superconducting state the constraints generate an asymmetric spectrum excitations between 
electrons and holes giving rise to an asymmetry tunneling density of states. 

Referee Comments: '"It is indeed refreshing to see an attempt at a completely novel 
route to some of these problems. The new approach presented in the manuscript has 
the potential of stimulating significant further developments by other researchers. I 
am looking forward for others to follow in the footsteps of the ideas presented in this 
paper" ' 
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I. INTRODUCTION 



The central problem in high Tc Superconductivity is to treat correctly the effects of 
strong electron-electron interactions. We consider the zero temperature region away from 
half fillings in the absence of the magnetic order. The physics in this regime is governed by 
the absence of double occupied sites. Based on experimental results we know that once the 
exchange interaction is added it will generate a superconducting ground state. 

For a lattice model the effects of interactions are described within a repulsive Hubbard U 
interaction. Due to the large on-site repulsion the double occupied state are prohibited. This 
means that we can project out from the electronic spectrum the double occupied states. As a 
result the ant i- commutation rules for the Fermionic operators are modified and calculations 
become difficult. 

A significant simplification takes place in one space dimension where the method of 
Bosonization shows that for any finite Hubbard U away from half fillings the physics is 
governed by two Luttinger Liquids (one for charge and the second one for spin). The limit 
of f/ ^ oo can not be considered in a microscopic formulation. This limit can be taken for 
the renormalized model within the Renormalization Group (R.G.) calculations. One obtains 
a line of Luttinger fixed points. 

n 

The wave function has been obtained from the Bethe ansatz [Ij in the limit U/t ^ oo 
(where t is the hopping constant). This solution shows that the spin configuration becomes 
degenerate at ?7 = oo and the wave function is a singlet ground state for all values of 
< U < oo. This means that the degeneracy at t/ = oo is removed by any infinitesimal 
perturbation t/U. A formal way for describing the wave function is to say that the ground 
state is annihilated by the singlet operator Q(x)\'^^ >= 0, Q{x) = ip'i{x)ilJi{x). The solution 
of this last equation shows that the wave function can be written as a product of a singlet 
state with another unknown state. 

For two space dimensions we do not have an exact solution as a function of the Hubbard 
U. Therefore we have to consider the projection of double occupancy which forces us to 
deal with the question of the modified commutators. A direct approach for dealing with 
the modified commutators is given by the Hubbard X operators. These operators are a 
mixture of Bosonic and Fermionic excitations. As expected, this leads to a complicated 
representation which is difficult to handle 21]. The conventional wisdom in higher dimensions 
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is the method of the slave particles [3| (slave Fermions, or slave Bosons). The slave particles 

for charge 
For space 



representation replaces the exclusion of double occupancy by two slave fieldsfone for charge 



and one for spin). These excitations are coupled by a U{1) gauge field [4, 
dimensions d > 2 the gauge field is in the confined phase, causing the slave particles to be 
strongly coupled. For some special conditions a deconfined phase might be possible s], ^ 
. The slave particles representation works in one dimension [ill] since the guage field is in 



the deconfined phase. In this phase, the excitations are described by solitons which carry 
fractional quantum numbers [lO|. An explicit solution based on the slave-boson method 
for the Hubbard U = oo case has been considered in the literature by [ll] within a path 
integral formulation. This formulation has been criticized in 12] which argued that the path 
integral measure for the slave particles is incorrect. The price we pay when we work with 
the slave particles is that the single particle spectrum is described by a pair of non physical 
excitationsl 

An alternative approach for dealing with the large U repulsion interaction is to use 

n n 

the Gutzwiller projection method used by p^, Using this method combined with a 
variational procedure, the authors [l^ have constructed a variational wave function for 
the strongly correlated superconductors. It has been pointed out ]15] that the exclusion 
of double occupancy is responsible for the strong asymmetry between the hole and the 
electronic excitations observed in the tunneling spectrum for the optimally doped BSCCO. 

One of the successful phenomenological theories used to explain a varieties of experiments 
is the marginal Fermi liquid theory la, llTI], yet the relation of this model to the microscopic 
theory is not clear. 

The recent quantum oscillations observed in the Shubnikov de Haas experiment might 
raise questions about the validity of different approaches and in particular it might test the 
validity of the projected wave function [l8|, [l9|. 

The purpose of the present paper is to introduce a new method for dealing with the 
problem of exclusion double occupancy. We propose to use Dirac's theory for First Class 
constraints [20]. The solution of the problem will be formulated in the language of Quantum 
constraints : for a hamiltonian H and a constraint operator Q{x) = ?/'|(a;)-?/'|(x) one has to 
find the many body state ]\E'° > which satisfies, if]\E'° >= £']\E'° > and is annihilated by the 
constraint operator Q{x)\'$^ >= 0. Since the constraint should be satisfied at any time the 
time derivative of the operator Q{x) requires that the commutator [Q{x),H] must vanish. 
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Away from half fillings we restrict the constraints to negative energies (holes excitations) 
and obtain a set of first class non-Abelian constraints Q'^^^x) a = 1,2,3 which obey 
Q'^~\x)\'^^ >= 0. In order to have a canonical theory in the presence of the constraints, we 



enlarge the Hilbert space by including new anti-commuting fields |21|, |22|. As a result when 
a hole is created at a time and destroyed at the time tf the evolution in the enlarged space 
is canonical. The Physical processes occur at times ti and t/ in the Physical Hilbert space. 
The projection into the physical Hilbert space is done by using proper boundary conditions 



for the anti-commuting fields [21|, |2^. Due to the non-commutativity of the constraint, the 
projection will generate a time dependent non-linear measure for the Lagrange multipliers. 
As a result, the physical evolution operator for the hamiltonian H with constraints is given 

where Ul=i1^K is a non — linear measure which is generated by projection into the 
Physical Hilbert space for the non-commuting constraints. The effective interaction for the 
holes sectors depend explicitly on the time interval T = — tf. 

The effective non local action obtained from the temporal projection is investigated with 
the help of the R.G. method. We find that the single hole excitation has a width which 
is linear in frequency and the scattering rate obeys ^ oc ujF{:^^) in agreement with the 



infrared data 



25( 1 . The linear frequency width is controlled by the holes density. When the 



hole density increases the region of the linear frequency width shrinks; for large densities it 
shrinks to zero and therefore a Fermi liquid behavior is obtained. The theory presented in 
this paper is applicable away not at half fillings for densities x > Xc at zero temperature where 
the magnetic order has been suppressed. Therefore we will not investigate the metal insulator 
transition. For a finite hole concentration the addition of a finite exchange interaction j ^ 
(in the t — j model) will give rise to a superconducting phase 26|] . An important question 
will be to understand the effect of exclusion of double occupancy on the superconducting 
state. We will show that by projecting out double occupancy, an asymmetry of tunneling 
density of states is observed. Using this theory we will explain the asymmetric tunneling 



density of states observed by [27 

The content of this paper is as 
method of First Class constraints 



bllowing: In chapter II we present the adaptation of the 



21 



22 



29|, l30|, l3l| to Condensed Matter Physics. In 



section III we present the model for correlated electron and show that the ground state can 
be obtained within the method of quantum constraints. The exclusion of the double occupied 
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states has been investigated in the hterature Q using the Gutzviller projection method. 
The solution of the equation Q(x)\'if^ >= is described by the Jastrow representation 

n n 

|13l . I14I ]. In chapter IV we show that the method of first class allows for additional non- 



Abelian constraints which give rise to a non-linear integration measure and replaces the delta 
function constraint used implicitly in refs. [isl. [l^ . Using the first class constraint method 
we compute the ground state for the high Tc superconductors. Chapter V is devoted to the 
computation of the integration measure. In chapter VI we introduce the canonical phase 
space action for the t — j model. In chapter VII we consider the explicit case where the 
exchange interaction is zero. In chapter VIII we introduce two Green's functions: G (the 
physical one) and D (the parametrical one). Due to the projection into the physical Hilbert 
space, the action which governs the single hole propagation is time dependent. The physical 
Green's function G is computed in terms of a parametric Green's function D. In chapter 
IX we perform the R.G. calculation for the effective action at the time interval T = ti — tf 
using the finite size scaling. Chapter X is devoted to the calculation of the parametric 
Green's function D. In chapter XI we present the calculation for the physical Green's 
function G. In particular we show that the relaxation rate is linear in frequency. We believe 
';hat this explains the experimental results observed for the optical conductivity given in 



251]. Finally in chapter XII we consider the effect of the projected Green's function on the 



superconductor ground state . We add a pairing interaction and compute the tunneling 
density of states for a superconductor using the projected Green's function. We show that 
the projection (incorporated into the calculation trough the self energy) give s rise to an 



asymmetric tunneling density of states, in agreement with the experiments 271. l28|. 

In chapter XIII we present our conclusion. We have included an Appendix where we 
discuss the effects of the secondary first class constraints generated by the commutators 
between the hamiltonian and the primary constraints. 



II. THE METHOD OF QUANTUM CONSTRAINTS 



The purpose of this chapter is to present an adaptation of the method of quantum con- 
straints to Condensed Matter Physics. 
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We have to find the ground state > under the conditions that a set of operators 
Qa{x), a = l,2....n are restricted to be zero! In Quantum Mechanics this means that one 
has to find the state > which is annihilated by the constraints Qa{x)\'^^ >= and is 
an eigenstate of the hamihonian >= >. 

There are two types of constraints : 

Second Class constraints are characterized by a non singular matrix 

[Qa{x), Qa'{x)] where the symbol [, ] represents the commutator for the Bosonic constraints. 
The matrix [Qa{x),Qa'{x)] can be inverted and therefore has a non vanishing determinant 
det[Qa{x),Qa'{x)] 7^ for a, a' = 1, ...m and m = |. Recently we have used this method 



to compute the persistent currents in coupled rings 23|] and study mesoscopic vortices in a 



two dimensional electron gas 24 1. 



When the determinant of the constraints vanishes we obtain First Class constraints 
which will be used to solve the problem of exclusion of double occupancy. According to 
Dirac [2^ one has to identify all the constraints which must be satisfied at any time. 

|Q«(f,t)|vl/0>=0 (1) 

From the Heisenberg equation of motion we obtain, -^Qai^x, t) = -r^lQa^x, t), H] where the 
new constraints are given by the difference between the commutator and a linear combination 
of the existing constraints, [Qa{x,t), H] — T^Q/3(a;)= Yln=it'^qf^{x). In this equation 
and stand for a set of matrix elements and qp{x), (3 = l,2..r represent the new 
(generated) secondary first class constraints. For the remaining part we will represent the 
two sets a = 1, 2..n (the primary first class constraints) and (3=1, 2..r (secondary first class 
constraints) by one set {Qi{x,t),Q2{x,t), ...Qn{x,t),qi{x),q2{x), ...qr{x)) =^ Qa{x,t) where 
a = 1, 2..{n + r) is the new index for the two sets. For the remaining part we will assume 
that Qa{x,t) are all the first class constraints (no new constraints are generated by higher 
order commutators of the hamiltonian with all the n + r constraints). 

Since the commutator of the constraints can be zero, the inverse of the commutator does 
not exists. As a result, a modification of the commutation rules as is done for Second Class 
constraints is not possible {3]. 

To overcome this difficulty one introduces new constraints $/3(x, t), (3 = 1, 2, ..(n + r), 

$^(x,t)|$°>=0 (2) 
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In order to obtain a canonical phase space [37j with less variables, we have to project 
out an even number of constraints (the constraints Qa{x) and their canonical conjugate one 
t)). This is achieved by demanding that the determinant of the commutator in the 
enlarged Hilbert space (with the additional unknown constraints $q,(s, t)) is not zero. 

Det[Q^, ^ (3) 
As a result we obtain an equivalent theory with a fewer independent degrees of freedom 



37(1 . Using this conditions we can compute the The Quantum Evolution Operator. 

a) The Quantum Evolution Operator for the unconstrained case is given by U[tf,ti] = 

The matrix elements of the evolution operator are computed according to the path integral 



method for Grassmann anti-commuting functions 38|]. Using the Grassmann Coherent states 
one introduces states |V'(t(3^) > ,< '?Ao-(a^)| which obey; iIj„{x)\%1)c,{x) >= ipaix)\ipaix) > and 
< ^ai^Mtix) =< tlj^{x)\^^{x) where {ip^{x)f = {ipa{x)f = 0. 

This allows us to formulate a field theory in the Schroedinger representation where the 
role of the coordinate is played by ipa{x) and the canonical conjugate momentum ipKx) is 
given by ipl{x) 



The matrix elements of the quantum evolution operator in the Grassmann space are given 
by < V',^/I^,^f_>=< i'\U[tf,tM >= til- 



Following 

U[^f,tf;ilJi,ti]: 



we obtain the path integral representation for the matrix element 



b) The Quantum Evolution Operator for the constrained system will be given in terms 
of evolution matrix elements Uphysl'ip f,tf;ipi,ti] in the Grassmann space. Due to the new 



constraints Faddeev 



391 ] has shown that a Super determinant is needed for the integration 



measure. Such a formalism has been used by 



40| 



A simpler method is to represent the physical evolution operator in terms of only the 
physical constrained Qa- The integration with respect to the unknown $^ constraints will 
modify the integration measure from dXa to a non — linear integration measure T>\a(T = 
ti — tf). This allows to represent the physical evolution operator in a form which is similar 



with to the result given by 



36| 



a=l 

Eq.(5) emerges from the canonical phase space formahsm which we will present in the 
remaining part of this chapter. 

We introduce the Lagrange multipliers, A"(x), a — 1, 2, ..(n + r) to enforce the exclusion 
of double occupancy. In addition we seek new constraints which are canonical conjugates 
to the original ones. The new constraints <l>^(x, t), (3 = 1,2, ..(n + r) are introduced with 
the help of new Lagrange multipliers T^a{x) ,ot = 1, 2, ..{n + r). The complete Hamiltonian 
which contains both type of constraints is given by the hamiltonian Ht'- 

(n+r) (n+r) 

Hr^H + jiY: A"(f )Q„(x) + E 7r„(f)$"(f)] d'x (6) 

a=l a=l 

It is convenient to replace the constraint by an equivalent constraint x"(x): 

^''{x,t)^-j^X"{x,t) + x''{x,t) (7) 

As result, the transformed Hamiltonian which contains the time derivative of the Lagrange 
multiplier A" (if) will be modified. 

Using the results given in equation (6) and (7) we obtain a new formulation of the 
constraint problem. At this point it is preferable to work with the canonical phase space 
momentum-coordinate action S. S — J dfix J dtL where L is the Lagrangian and h{x) is 
the hamiltonian density for the hamiltonian H. 

(n+r) (n+r) 
in^iix,t)dtMx,t)+ 7la{x,t)dtX''{x,t)-h{x,t)- Y {T^a{x,t)x''{x,t)+r{x,t)Qa{x,t)) 

(T=T,i a=l ot=l 

(8) 

The term 7ra(x, i)5tA"(x, allows to identify the canonical conjugate momentum of A"(x) 
with 7ra(x) which obeys the commutation rules, [\°'{x),T:p{y)\ — ih5a,p5{x — y). 

It is important to point out that equation (8) can be understood as a starting point 
for our theory where the Lagrange multiplier \"-{x) enforces the constraints Qa{x). The 
Quantum nature of the Lagrange multipliers is enforced by demanding the existence of a 
canonical conjugate variable Tia{x)- Physically the canonical conjugate momentum must be 
enforced to be zero. This is done by introducing a new Lagrange multiplier x"(x). As a 
result our theory will have the set of constraint equations: 
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7r„(f)|^° >= and Qa{x)\^° >= 0. 

Which will be enforced by the two sets of Lagrange multipliers: ^"^^ A"(x). 

The new Lagrange multipliers x'^i^), oi = 1,2, ..(n + r) have not yet been specified. 



VILKOVISKY 



29| has proved a theorem which shows that the role of the new Lagrange 



multipliers x'^jx) is eq uivalent to the gauge fixing in Quantum Electrodynamics. The proof 



of the theorem 



29 



30( 1 is based on the followings steps: 

1. The action given in eq.(8) can be represented in an equivalent way using new fermionic 
fields. One introduces a pair of anti-commuting fields for each one of the constraints 
7r„(f)|^° >= and >= 0. 

2. The Lagrange multiplier x'^{x) plays the role of gauge fixing condition. The replace- 
ment of the action in equation (8) with a new action written in terms of the new Fermionic 
fields show explicitly that the expectation values for any physical observable is invariant un- 
der a change of the Lagrange multipliers x'^ix) — ^ x'^ix) + ^{x°'{x))- Therefore the physical 
results are independent of the particular choice of the fields x°'{x)- 

An exact mapping of the action given in equation (8) to an equivalent action 



29|, 



3Q|. 



where the constraints are replaced by new Fermionic fields is possible [22 
As a result one obtains an enlarged Hilbert space of anti-commuting fields. 

We consider the case where the constraints Qa{x) obey the following relations: 

[Qo.,Qp\ = Flp{Q,) (9) 

Where F^^piQ-y) is a nonlinear function of the constraint such that for = we have 
the result Fl^piQ-f = 0) = for any a, /5 and 7. 

For this case the VILKOVISKY [29] theorem allows the exact mapping of the action 
in eq.(8) to the new action defined in terms of the new anti-commuting fields. 

The mapping is done according to the following steps |22| : 

a) For each constraint field Qa{x) one introduces a pair of anti- commuting real fermions 

C°(x) and ba{x) 

[C"(x),6^(x')]+ = So,,pS{x~x') (10) 
Such that bfs{x') acts as an annihilation operator. 
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b) Similarly the canonical momentum constraints 7rQ,(x) is replaced by the pair of anti- 
commuting real fermions e°'{x) and fa{x)- 

[e"(x),/^(x')]+ = 5a,/35(f-x') (11) 
c) The physical Hilbert space is extended to an enlarged Hilbert space. Therefore 

the many body state is replaced by a new state |\1' > which is build from a Fermionic 
subspace needed to enforce the constraints. The wave function < tlj„{x)\'$^ > is replaced by 
a wave function in the enlarged Hilbert space 

< V'.(x),C"^(5');vr„(f),e"(x)|* > (12) 

Where < ipa-{x),C°'{x);Ha{x),e"'{x)\ is the coherent state representation in the enlarged 
Hilbert space. 



d) According to the theorem (see pages 247, 322-324 in (2lJ) the physical wave function 
< 4'„{x)\'^^ > is obtained by the projection < ip„{x),C°'{x) = 0;TTa{x) = 0, e"(x) = 0|\1/ >. 

The physical evolution operator matrix elements Uphys['ipf,tf;ipi,ti] are obtained 
once we impose the temporal boundary conditions on the auxiliary fields in the 
enlarged Hilbert space. We use the following temporal boundary conditions: 

C°(f,t,) =C"(f,t/) = (13) 

e"(f,ti) = e"(f,t/) = (14) 
For the momentum 7ia{x,t) we use the following boundary conditions: 

Tlaix, ti) = Tlaix, tf) = (15) 

e) In the extended Hilbert space the wave function > is replaced by |\E' >. 

In the enlarged Hilbert space we replace the constraints Qa{x) and liaix) by a 
new constraint Q. 

The constraint Q operator has to obey 

n(x)|^>=0 (16) 
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f) In the extended Hilbert space the hamiltonian H is replaced by He- The operator Q, 
obeys the extended Heisenberg equation of motion: 



in^Vt\^ >=[Vt,H^]\^ >=Q (17) 
The extended constrained operator Q must be NILPOTENT 

Q2 = [n,f}]+ = o (18) 

g) In order to satisfy the NILPOTENCY condition the operator Q, is written as a sum of 
two parts 

O = Oo + ^NL (19) 

Where is given by 

(n+r) (n+r) 
a=l a=l 

Due to the non commutativity of the constraints [Qa,Qp\ — F2,i3iQi) con- 
dition = requires that the constraint operator should have a non linear pcirt 
given by I^tvl- The most general form of the non linear part VLml is given by: 

^nl{x) = EtY^ c-m d^x' EStf E^:;? .. / d^x^-^^) E£t? ^ 

• (f ^)..C"(-'^) (f (-+'-))Mf-(^-;^6^, (^')-^„,., (^^"+^^)] (21) 
The matrix Mf?, v^"+'^ elements are determined by the condition VL^ = 0. 

h) The hamiltonian He represents the extention of H for the extended Hilbert space: 

/{n+r) (n+r) 
d^x E E C'i^KixMx) + ... (22) 
a=l ^=1 

The parameters h^{x) are determined by the equation [Hf., il] = 0. 

k) The condition >= guarantees that flOphysl"^ i^TOphysW"^ >■ As a result 
any physical operator Ophys is mapped by the operator Q, to another physical state, therefore 
we have [O,, Ophys] — 0. This implies that if |^ > is a physical state, then QOphys]^ >— 

[n,Ophys]\^ >= 0. 

11 



This condition shows that any modified physical operator Op^^^ which is related to the 
original physical operator Ophys trough the transformation Opf^y^ = Ophys + [^,W]+ has 
the same matrix elements as the original operator. The symbol [, ]+ stands for the anti- 
commutator and ly is a new fermionic operator defined in terms of the Lagrange multipliers 
A°(x) and x'^i^) the Fermionic fields. 

(n+r) 

W = d'^x [&a(x)A°(f) + ea{x)x'"{x)] (23) 

The proof that the operator O'pf^y^ = Ophys + [^5^^]+ and the operator Ophys have the 
same matrix elements in the extended Hilbert space follows from the identity [Q, [Q, W]] = 
[(]2^ W] = (the nilpotency of fi). 

m) This means that the matrix elements of the hamiltonian H effective are the same as for 
the hamiltonian H^.. 

Heffeetive = + [Q, W]+ = + (24) 



This shows that 6q~^ acts as an exterior derivative [2l|, 6'q'A = [Q,A]^ where A is an 



arbitrary operator. 



n) Following the theorem [29|] we have the freedom to choose any bosonic fields Iii 
particular the path integral is independent on the choice of the unknown Lagrange multiplier 
X°(5'). The path integral is invariant under the transformation x'^i^) ~^ ix^ix))' = X^i^) + 
5(x"(x)). 



As a result of this theorem one can show 



2l| that the quantum expectation value of any 



operator A{t) which commutes with the constraint [A, Q] = is independent of the choice 

< Mi) >x"=< Mt) >x-+Hx-m (25) 

The method presented in this section will be used to solve the problem of exclusion of 
double occupancy in the next chapters. 
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III. THE t - J MODEL FOR CORRELATED ELECTRONS 



In this section we will present the model for the high Tc Superconductors. We will show 
that the ground state can be computed using the method of First Class constraints 2o| . 

For the case that the hopping parameter t and the one site repulsion U obey the condition 
(t/U) < 1, the double occupation states are projected out and one obtains an effective 
hamiltonian H =^ PHqP + ■P^o(i--P)(i-^)goJ' _|_ where P = 1 — rijnj is the projection 
operator. The effective model is given by: 

H = -tJ2 J2 ^t{x)M^ +a) + h.c. + 6H = H0 + 6H (26) 

X,d (7=1,1 

where x represents the lattice points and a runs over to the nearest neighbor sites. This 
is the t — J model where Ho is the hopping hamiltonian which acts on a restricted Hilbert 
space where double occupancy is excluded, while 6H represents the exchange hamiltonian 
controlled by J oc ^. 

The Many Body ground state is given by |\1''^ >. The exclusion of double occupancy on 
each lattice point x is imposed by the constraint condition which determines the ground 
state >: 

V^Tm(x)|vi>°>=0 (27) 

We define the constraint field Q{x) = ■i/ji{x)iIji{x). To find the ground state |\E'° > of the 
hamiltonian Hq + 6H which obeys the constraints we have to satisfy the following equations: 

i7|^°>=E|^°> and g(f)|^° >= (28) 

Once the ground state |\E''' > is found, the excitations spectrum is obtained by applying 
the creation operators on this ground state such that no double occupied excited states are 
created. 

It is easy to see that the solution to eq. (28) can be written in the form: 

>= n,^[l - Q+{x)Q{x)]\^ > (29) 
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Where > is a state which must be determined! The state > belongs to the 
13] and is often determined by variational methods {l4]. Two possible 



Gutzwiller class 



choices for |\E'° > have been considered: the BCS wave function has been used to compute 
the RVB state and the Fermi Liquid ground state \F.S >= U^'^^'ip^ {K)tp~l^ {K)\0 > has 
been introduced to compute the strongly correlated metallic state. 

In the present paper we will not perform a variational calculation, instead we will compute 
the ground state wave function by determining the additional conditions which the state 
|\E'° > has to satisfy. We will show that using the Fermi surface as a the unperturbed ground 
state we find two additional constraints, the pair creation constraint Q^{x) and the hole 
number constraint Qsix) = l—n^{x)—ni{x). The set of the three non-commuting constraints 
replace the delta function measure (which results from the condition Q{x)\'$^ >= ) by a 
non-linear integration measure . It is this measure which generates an effective temporal 
interaction with a time dependent coupling constant. 

These results will be derived in the next chapter using the method of first class constraints 
introduced in chapter //. 



IV. THE APPLICATION OF FIRST CLASS CONSTRAINTS TO THE PROBLEM 
OF EXCLUSSION OF DOUBLE OCCUPANCY 



In this chapter we will apply the general theory for First Class constraints introduced in 
chapter II to solve the model presented in chapter III. The eigenfunction |\I'° > must obey 
Q{x)\'^^ >= 0. We will show that away from half fillings we can use a set of first class 
constraints defined only for negative energies Q^\x) where a = 1,2,3. The commutator 
[Q^^\x),H] generates secondary constraints qlr\x), r = 1,2,3 which will be neglected for 
describing the low energy physics. The justification for this approximation is given in the 
appendix of this paper where we show that the effective action generated by the secondary 
constraints are irrelevant according to the R.G. analysis at low energies. 

Due to the constraints, the non interacting Fermi energy will be shifted by Sfip to a. new 
value. The metallic behavior will be characterized by the vanishing of the renormalized 
chemical potential shift = [431]. 
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At half fillings we have two additional constraints: Q^{x) and the hole number opera- 
tor Qsi^x) = 1 — [ip]^{x)ilj^{x) + ip]^{x)'ilJi{x)]. At half fillings the three constraints satisfy: 
>= 0, >= and (53(f) 1^^° >= 0. Away from half fillings we have 

only one first class constraints Q{x), the other two constraints are neither first class nor 
second class. This difficulty can be resolved by modifying the constraints. Away from half 
fillings we will restrict the constraints only to negative energies /holes excitations Q^~\x) 
where a = 1,2, 3. We introduce the definitions: 

MS) = i^l'Ks) + i^l-\xy,Ms) = ijt^+\x) + ^t^-\x) (30) 

The notation (+) represents the particles excitations for positive energies and (— ) de- 
scribes the holes excitations for negative energies, both measured with respect to the renor- 
malized Fermi energy. 

tp^^\x) and ^}}~\x) act as a destruction operators with respect to the ground state 
> and obey V^(+)(f)|vl>o ^t(-)(f)|\[/0 >= Q. From [41] we learn that the field M^x) 
is build from the Fourier momentum components K < Kp, where Kp corresponds to the 
non-interacting Fermi momentum. Similarly il)^\x) is build from the Fourier momentum 
components K > Kp- The effective renormalization of the Fermi surface caused by the 
constraints is included into the chemical potential shift 6fip, which will be computed within 
the R.G. calculations. 

The renormalized ground state > has no particles above Kp and no holes below Kp. 
(At this step we do not make any assumption about the nature of the discontinuity of the 
Fermi Surface (the occupation number at Kp). 

Acting with the constraint operator Q{x) on the ground state > we observe: 

Q(f)|^° >= {^\'^\x) + ^\~\x)){^['^\x) + ^j"^(f))|^° >= Vj"^(f))^{")(f)|^° >= (31) 

This equation shows that the constraints can be restricted to the holes type excitations 
Q*^~-*(x) = V'j ''(f))V'j \S)- For positive energies the constraint is automatically satisfied, 
therefore the constraint Q{x) is restricted to Q^^\x). 

Using the holes representation, ijjl~\x) and ipl^~\x) we construct the new representations 
for the constraints: 
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Q^-\x) = ip\~\x)il;'f\x) (32) 

Qt(-)(f) = [^^{-)(f)^^{-)(f)]t (33) 

g^-)(f) = 1 - [ilj\^-\x)i/j\~\x) + i/jf'\x)ijjl\x)] (34) 

It is convenient to work with real constraints. We introduce two real constraints, Qi"* (x) 
and Q2 \x) using the constraint Q^~\x) and Q'^''~\x). 

QV{x) = -^iQ^-\x) + Q^^~\x) (35) 

Q^f\x) = ^(g(-H^) - Q^^-H^)) (36) 



In addition to this two real constraints we have the third real constraint (x) given by 
equation (34). The set of the hole type constraints (5^~-*(a;)|\l/° > ,« = 1,2,3 represent the 
Non-Abelian First Class constraints for our problem. The commutators of the constraints 
obey: 

[Qlr\s),Q'^i\s)]-ira,pQVi^) (37) 

The Non-Abelian First Class constraints are characterized by the structure constants 

/J,/3' f2,/3 = ~f}.a = 1 for a 7^ /? 7^ 7 and zero otherwise. 

The commutator of the constraints with the hamiltonian H satisfies the equation: 
[Q^a\^)i ^\ ~ '^p(^)Qp~\^)i where Tp(x) stands for a set of local matrix operators. The 

symbol means that the derivatives of the constraints operators have been neglected. This 

approximation can be seen by computing the commutator of the constraint with the kinetic 

energy: 

['4ja=l{x)ilJa=^{x),Ho] = tEa[^a=i{x)ipa=^{x + a) + iIj„=i{{x + a))ip^=i;{x)] 
The new secondary constraints are obtained by taking the difference of the commutator 
^[Q{x), Hq] with the primary first class constraints Q{x): 

q{x) = Ea[^<T=i(^)^/'<T=T(^ + a) + iJa=i{{x + a))V^<,=|(f)] - Q{x) 

q+{x) = T,a[fPa=l{x)lpa=^{x + o) + V'a=i((:? + a))V'(T=T (^)]'*" " Q'^ (x) 

and 
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qsix) = T,abPa=l{x)tpa=l{x + a) + 1pa=^{{x + a))7jja=f{x)] - [Qsix) - 1] 

Using the hnear transformation given by equations (35) and (35) we obtain the new set 
of secondary first class real constraints: ql~\x), r = 1,2,3. 

The secondary first class constraints modify the Lagrangian L hj 6L : 
SL = El=i Mx, t)dty{x, t) - j:l=i{Mx, t)r{x, t) + X'{x, t)qi-\x, t)) 
where A''(x, t) are the Lagrange multiphers which enforce the secondary constraints and 
7tr{x,t) are the canonical momentum conjugated to the new Lagrange multipliers y'{x,t). 
The integration over the new Lagrange multipliers will generate products of q^~\x) pairs. 
Such operators will induce effective interactions which have the dimensions of the square 
of the kinetic energy with two time integrations. According to the analysis given in the 
Appendix the engineering dimensions for such operators is —1 and therefore they arc strongly 
irrelevant at low energies ( in comparison to the engineering dimensions 1 for the effective 
interaction due to the primary constraints) . Therefore we will ignore the Lagrangian induced 
by the secondary constraints or/and higher order constraints generated by the commutators 
[ql.~\x), Hq], r = 1, 2, 3 (such operators contain higher order derivatives in comparison with 
the constraints ql.~\x)). 

The canonical phase space action for the holes constraints which replaces eq. (8) (without 
the generated constraints described by SL) is given by: 

L^ ^hi^l{x,t)^tMx,t)+j2^a{x,t)^tX^{x,t)-h^^ 

(38) 

where 7ra{x) is the canonical momentum conjugate to the Lagrange multiplier A"(x). 
This action has two sets of constraints (5^~^(x)|^° >= and 7ra{x)\^'^ >= 0. The Lagrange 
multipliers are: A°(x) and the unknown one x'^ix)- Using the general theory for First Class 
constraints presented in section II, we will express the action in eq. (38) using new anti- 
commuting fields. The anti-commuting fields C°'{x) and ba{x) are used to replace the holes 
constraints Q^~\x^t) and the anti-commuting fields e°'{x) and fa{x) replace the momenta 
constraints naix)- The unknown Lagrange multiplies act as gauge fixing conditions. 

We replace the state > by the state |^ > defined for the extended Hilbert with the 
extended constraint operator Vt for the extended Hilbert space (see eq.(19)), = fig + ^nl 
where Qatl is chosen such that the condition = is satisfied. 
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The operator Qq (see eq.(20)) is given in terms of the physical constraints \x) and 



On 



d'^X 



Ec"(f)g(-)(f) + Er 



[X n^ix 



.a=l 



a=l 



(39) 



The Nonhnear operator Qnl (see eq.(21))is obtained using the structure constants given 
in eq.(37). 



^iVL(x) = EC"( 



X) 



a=l 



3 3 



d'^' E E -irCHx)fl,b,{x')] 

/3=l7=l 



(40) 



According to the theorem |29|] the path integral is invariant under the change of the 
Lagrange multiphers x'^ix)- For the continuation of this article we will choose x'^ix) = 
and obtain from eq.(23) the Fermionic operator W. 



3 

W= d'^xJ2 baix)\ 

•' a=l 



[X 



(41) 



We compute the anti-commutator and find: 



[w, nu = E +thj: h,{x)nx) + lc^ix)nx, t)fi,b,ix, t) 

a=l a=l a=l 13=1 "1=1 

(42) 



This result allows to obtain the effective hamiltonian in the enlarged space: 

Heffective = He + [fi, W]+ whcrc He is choseu such that [He, fi] = 0. 

In this case the difference between He and H is given by operators which contain higher 
order derivatives and therefore are irrelevant for low energy Physics. 

Using the effective hamiltonian HeffecUve = H + \p.,W]^ we obtain the action as a 
function of the initial and final times ti and t/, T = ti — tf . 



s[t„tf] = f' dt[ E Hh)^i{x, t)dtMx, ^) + E ^-(^' t)dtX''{x, t) 

<T=l,t a=l 



Y,dtC''{x,t)be,{x,t) + Y,dtee,{x,t)rix,t) - 



a=l 



a=l 



[h{x, t) + E A°(^, t)Ql-^ (f , t)+zhY^ 6„(x, t)r (f , t) + 



a=l 

3 3 3 1 



a=l 



(43) 



a=l p=l 7=1 
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For commuting constraints, the last term in eq.(43) vanishes (the anti-commuting fields 
are decoupled). For such a situation the constraints in eq.(43) can be treated by delta func- 
tion constraints. Once the non-commutativity of the constraints is considered one generate 
a non-linear measures for the Lagrange multipliers (due to the last term in eq.(43). 

The full Non-Abelian action given in equation (43) describes a canonical evolution in 
an enlarged Hilbert space. The physical observables are obtained only after we project 
the extended state into the physical Hilbert space. This is achieved by imposing temporal 
boundary conditions on the non physical anti-commuting coordinates C"(x, t) ,ea{x,t) and 
the conjugate momenta Haix, t) ,a = 1, 2, 3. 

C°(x, ti) = C°'{x,tf) = ea{x, ti) = ea{x,tf) = na{x, ti) = 7la{x,tf) = (44) 

The projection into the physical space is achieved by the use of initial time ti and final 



time tf boundary conditions 211 ]. 



Uphys[iJfitfy,^^J^iU)] = U[^f{tf),TTaix,tf) = O^Cix^tf) = 0,eaix,U) = 0;^i{U)] (45) 
The physical evolution matrix elements are given by: 



UpHys[i^f{tfy,ij,{U)] = J D(7/;^,V^.)D(C",6jD(e„,r)D(A-,7r„)exp[^(V^^V^;+^,V^,+^^[t„t/]] 

(46) 

Where S[ti,tf] is the action in equation (43) . D{ip^,'(p^) represents the Grassmann 
measure, D{C",ba) and D{ea, f") represents the integration measure j2l| for the fictitious 
non-commuting degrees of freedom. D{X'^,TTa) is the phase space integration over the La- 
grange multipliers and the canonical conjugated variables. The action S[ti,tf] can be written 
as a sum of the physical action S[ti,tf] and a non — physical action S^^°'^'^[ti,tf] given in 
terms of the non-commuting degrees of freedom C^fia and Ca , 

S[ti,tf]=S[t„tf] + S^''°'%,tf] (47) 

For fixed values of ti and we integrate over the non— physical non-commuting degrees of 
freedom C^fia and Cq , We obtain an effective action S''^-^'^ [A"; ti,tf] as a function of the 
time interval and the Lagrange multipliers. This action is obtained once the non-physical 
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fermions fields have been integrated out and the boundary conditions for the conjugate 
momenta TTaiU) — Tra{tf) — have been used. 

^.5e//.[A«;*.,,] ^ J i^(C",6„)D(e„,r)ei^^"^*^>*/lU„(,,)=.„(,,)=o (48) 

This allows to represent the matrix elements of the physical evolution operator by the 
formula: 

U,nys\^f(tfyMti)] = / L'(V^,,V'.)i?(A'^)e^^^''-[^"^*-*^le[i(^/^/+^^^')+t^[*-*/]] = 

J D{^^, ^a)I^(A"; U - t^)e[5(^/^/+^*^*)+i^[*-*/ll (49) 

where T>\{x] U — tf) is a time dependent non — linear measure obtained by integrating 
out the non-physical anti-commuting fields which describe the dynamics of the canonical 
conjugate fields. The measure is defined trough the path integral in eqs. (48)-(49). As a 
result any physical Green's function will depend on a coupling constant which is a function 
of the time interval T — ti — tf (generated by the projection of the non physical Fermionic 
degrees of freedom at the initial and final time). 



V. COMPUTATION OF THE INTEGRATION MEASURE P(A"; ti - tf) 



The physical evolution operator given in equation (49) depends on the integration measure 
T>{X°'; ti — tf). This measure is computed from the action S^'^°^*[ti,tf] and is defined by the 
equations (43), (47)-49. The exphcit form of S<^^°'%,tf] is given by: 

'-S<^^"^%, tf] = [f: 9,C"(x, t)h^{x, t) + j2 dteaix, t)r{x, t)-j2 K{x, t)rix, t) + 

n Jti Q,^^ ^^-^ ^^-^ 

E E E \cP{^)nx. t)f2,,b,{x, t)] dt (50) 

a=l l3=l 7=1 

The integration over the non-physical fermionic field is done using the variation of the 
action 5[S^^"^^[[ti, tf]] = 0. We obtain the relation dt{ea{x, t)) = ba{x, t). As a result we find 
that the integration measure is given by V{\'';ti,tf) = D(A°)e^^'^^[^";*-*/l. The effective 
action e^*^"^^ '^"'^''^J"! is determined by the Grasmann integration: 
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ICdl-/-^ - J-/oV««)l7)) (51) 



with / = I{a,f3) being the identity matrix. Following [45|, |46| we evaluate the two deter- 
minants in eq. (51). The determinants are a function of the time intervals T = ti — tj. 
For convergence reasons at T — oo we rotated the time contour [42]. We define a new 
time interval f = Te'^. The effective action S^ff-[X^]tie'^,tfe'^] = S^^f-[ip{x,f);tie'\tfe'^] 
is given by, 

S'ffM^,f)-U,tj] =A' f d'xLog[7r^^!^^^^i^] (52) 



ip{x, T) 

where A is the ultraviolet cut-off and d is the space dimension. The integration variable 
A°(x) is replaced by ip°'{x,T) given by (p°'{x,T) = j^^i^M with the amplitude (p{x,T) = 



f v/(^^)2 + (Al^)2 + (A(fM)2. Eq.(52) represents the new integration measure for the 
Lagrange multipliers. 



VI. THE EFFECTIVE MODEL FOR EXCLUSION OF DOUBLE OCCUPANCY 
IN TWO DIMENSIONS 



The action in equation (43) can be written in a simplified form once we replace the 
measure dX°'{x,t) with the non-linear measure d\'^{x)e v{3,t) 

S[t,,tf] = /// dt[Y.„^^,^{-th)i;l{x,t)dti;,{x,t) - [h{x,t) + Yi=i\''{x,t)Q^-\xM 

The new measure d\'^{x)e v(s,t) -^as obtained as a result of the non-commuting 
constraints Q^~^ (x, t) . The change in the measure can be rewritten explicitly in terms of the 
effective interaction given in equation (52). 
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VII. THE EFFECTIVE MODEL IN THE ABSENCE OF THE EXCHANGE IN- 
TERACTION 



In this section we will present the effective model for a free electron system which obeys 
the exclusion of double occupancy. 

We will consider the case where the exchange term 5H = 0. Therefore wc will replace 
h{x,t) ho{x,t) (in equation (43)). We believe that this model describes the situation at 
zero temperatures away from half fillings above certain holes concentration Xc (where the 
magnetic order is absent) . When the exchange is zero we will not be able to investigate the 
metal insulator transition, which takes place close to the half filled case x — > 0. We will 
consider the case where the holes density x obeys x > x^. Experimentally we know that 
at T=0 such a window exists between the magnetic ordered state and the appearance of 
superconductivity. In order to study the nature of the metallic state at zero temperature and 
holes densities x > Xc it is reasonable to suppress the superconducting order by restricting 
the exchange interaction to zero. We will show that at low holes concentrations Xc < x « 1 
the imaginary part of the self energy is proportional to cu. For large holes densities x ^ 1, a 
crossover transition to a Fermi liquid with the imaginary part of the self energy proportional 
to LV^ is obtained. 

In order to perform a R.G. study we use the free Fermi liquid representation in two 
dimensions. 

S'[ti: tf] = Jo ¥ /I dt I %J[^^ E.=T,i[^i(^> - ne)R,{t, e; s) + 4(i, e; s){indt + 

he)L^{t, e; s) 

where R(r{t,e;s) and Lcr{t,e;s) are the right and left chiral fermions in the channel s. 
The polar angle s on the Fermi surface is restricted to the region [0 — tt]. The momentum 
excitation normal to the Fermi Surface is given by where vp is the Fermi velocity. J[e, s] — 
^l^^^g)! is the Jacobian transformation from the K^; and Ky coordinates to the energy e and 
polar angle s. We will rescale the Fermionic field by Jo ^J['^, s]. As a result, the Fermi 
surface action will depend only on the dimensionless Jacobian J[e, s] — r- ^ '^^'f, — y- 

The effect of the constraints are supposed to shift the position of free Fermi surface given 
by the Fermi vector Kp[s). This effect will be taken in consideration by a finite shift of the 
chemical potential 5^f{s). As a result, the non interacting action S'\ti,tf\ will be replaced 
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by S'''[tj,t/]. The complete action for our model as obtained from the previous chapter, 
including the shifted Fermi Surface, is: 



S[U,tf] = S%,tf] + S'-''"'%,tf] + S'"-[ip''{x,f);U,tf] (53) 

where S^ff-[X'';ti,tf] is given in equation (52) and represents the modification of the 
integration measure for the Lagrange multipliers A"(a;). S^~^°''^[ti,tf] which originally was 
given by Y.a=i A"(x, t)Q^^\x, t) is written in terms of the new variables ^Pa{x, T), a = 1, 2, 3: 

S^-'^''^%„tf] = t dt [ d'x[j: ^^^§^g(-)(f,t)] (54) 
S%Sf] = T- r dt I ^J[e,s] [Ri{t,e;s){zhdt-he)K{t,e,s) 

Jo 71 Jti J ZTT t 

+Ll{t, e; s){ihdt + he)L„{t, e; s) + 

Rl{t, e; s)6fiF{s)Ra{t, e; s) + Ll(t, e; s)SfiF{s)L^{t, e; s)] (55) 

where S^[ti,tf] represents the free fermion action and 6fiF{s) is the shift in the chemical 
potential induced by constraints. The metallic phase will be identified by the vanishing of 
the renormalized chemical potential shift [5^F{sy\R = 43 1. 

The action in eq.(53) represents the complete solution for the problem of exclusion of 
double occupancy for the tight binding fermion model. It is important to mention that the 
action describes an effective interaction for the holes which depends explicitly on the time 
interval T = ti — tf. We have a situation where the Lagrange multipliers A" (a;) can be 
treated as random variables, similar to the situation for annealed disorder in Statistical 
Mechanics. 

Next we compute an effective interaction which is induced by the constraints. We in- 
tegrate the Lagrange multipliers <y9"(x, T). We perform the integration for times T > 
where A is the momentum cut-off and vf is the Fermi velocity, which is maximum at half 
fillings and decreases with the doping. The product f^A represents the electronic bandwidth 
in frequency units. Keeping only terms which are second order in <^"(x, T) in the action, 
Loqlir allows us to integrate out the Lagrange multipliers and obtain the effective 
interaction for times larger than T > : 

a) For electrons no interaction is generated and the physics is given by: 
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Jo TT Jti J ZTT 37, 



+4(+)(t, e; - he)Li+\t, e; s)] (56) 



b) For holes we find the following action is generated at long times T > 



S^ff-^'^^%,ff] = T- fj dt I ^J[e,s] [Ri^-\t,e;s)(ihdt + he)R^-\t,e;s) 

Jo 71 Jti J ZTT . 

+Ll(-\t,e;s){indt + ne)L';;-\t,e;s) 

+Rl^-\t, e; s)SMs)Ri~\t, e; s) + L'^y\t, e; s)5iJiF{s)L^{t, e; s)] 
+ [ — I — [J dti [J dt2 n 7^^'^[ei,si]<^[e2,S2] 

JO TT JO TT Jti Jti (ZTTj'' 

- S2; f)5{-ei + 62 + 63 - 64) 

[[Rf-^ (ti, -ei; si)Lf (ti, €2; (ts, eg; S2)i?["^ (ts, -64; ^2) 

+Lf-)(ti, ei; Si)R\'^'\h, -62] Si)L\-\t2, £3; S2)i?(-)^(t2, -£4; S2)] 
-ei; si)L|("^(ii, 62; si)i?{~^(i2, -63; S2)4"^(i2, €4; S2) + 
Lf~\ti, ei; -62; si)4-^(t2, 63; S2)i?{"^(i2, -64; S2)]] 

^ 5o'°^^[^ , 0] + S!il-''''%, tf] (57) 

For a spherical Fermi surface we have J[e, s] = 1 = e° and for all other cases we have 
J[e, s] ~ where < 2; < 1, where A is the momentum and energy E = vpA cut-off's. 
This allows to represent the coupling constant by g{s-T) = g[s;T)vFE. We will use the 
dimensionless time t — = ^ and define a real and imaginary coupling constant: 

g{s; T) = —i{u — i/S) where the u = u{t) — 3( y^^^ )^ and A is given by A = uSin{5) where 
5^0. The effective action Slll~^"^%Jf\ in eq.(57) resemble the hole- hole interaction 
for Superconductivity at different times. This form is obtained after we have replaced the 
constraints operators Q[^\x,t) by the chiral fermions Rfj{t,e; s) and L„{t,e; s). 

Next we consider the situation away from half fillings, without the umklapp terms. There- 
fore we will not attempt to describe the Metal Insulator transition for which the presence of 
the exchange and umklapp interactions are important. The action obtained is valid above 
a certain critical hole concentration x > Xc- The theory derived depends on the bandwidth. 
The bandwidth decreases with the increases of the hole concentrations. Experimentally it 
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is observed that for large hole concentrations the Fermi Liquid behavior is observed again. 
At this stage we can not prove that we have a transition to a Fermi Liquid. We can show 
that with increasing holes concentrations, the frequency region for which the self energy is 
linear in uj shrinks to zero. 



VIII. THE SINGLE HOLE GREEN'S FUNCTION 



The single particle excitations are the same as for the free electron case. The only change 
is that the holes type excitations are governed by the action given in equation (57). 
The definition of the retarded Green's function for holes is given by: 

G^-\x,a,tf-x',a\U] = id[U-tf] < ^ViP(x', t^^l'^^, > (58) 

where d\ti — t/] is the step function which is 1 for time intervals ti — tf > and oth- 
erwise. An alternative representation for the Green's function is possible if we use the path 
integral representation 37|]. In particular the matrix elements of the evolution operator 
UphysH' f,tf;ipi,ti] in the Grassmann space [371] allows us to compute the Green's function. 
The relation between the path integral representation and the direct computation of the 
Green's function G^''^x,a,tf]x',a',ti\ using the ground state |\E''' > has been shown in 
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42| . In our case, the effect of the time dependent coupling constant can be investigated 



within the time dependent path integral given in the Grassmann space, as below: 



\ilji{x',a',ti),tf,ilJi{x,a,tf),tf] = i <ilJi{x',a',ti),ti\-^l'^, \x,tf)ilji \x,tf)\iljf{x, a,tf),tf > 

= I I D(^,,V^.)V^,(x',a',t.)V^;(f,cr,t^)e^^''"''°"[*-*/l (59) 



Following 
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42| we project the j{x\ a',tf),tf] 4'i{x, a, tj), tj] into the ground state 

> and obtain the physical Green's function G^~''[x',(T',tf;x,(T,ti] defined in eq. (58) 
which is represented in terms of the matrix elements in the Grassmann space. 

<^,{x\a',ti) = 0,t,\ijl')~\x,tM-\x,tf)\^f{x,a,tf) =0,tf > 



[x', a ,tf; X, a, ti] 



i < i)i{x\a',ti) = 0,ti|?/'/(f, a,tf) = 0,t/ > 



< F.S.\^l\-\x,tM-K^.tf)T[e^ 



i qeff-hole 



M]\F.S. > 



< F.sme^'tJ:ir'°'%^mF.s > 



(60) 



25 



where \F.S >= H^'Lq {K)\0 > represents the Fermi Surface and T [...] represents 

the time order. 

The Physical Green's function G^~\T, K;g(T),T, A] will be computed using the finite 
size action S^-f^~^°''^[ti,tf;g(T),T,A]. The finite size effect is introduced by the duration 
of the hole excitations ti — tf = T. The projection at times ti and tf generates a time 
dependent action. Due to the fact that the coupling constant g{T) is a function of the time 
duration for the hole excitations, it is advantageous to introduce a parametric Green 's 
function D^^^[x' , a' ,Tf; x, a,Ti] with the coupling constant g(T) and the parametric time 
interval t = Ti — tj where, tf < Tf < Ti < t^. 



n(-)r^, / 1 _ <i^ii^','^',U) =0,U\ijli \x,Ti)tljl \x,Tf)\ijf{x,a,tf) =0,tf> 

i < ipi{x', cr', ti) = 0, ti\ipf{x, cr,tf) = 0, tf > 
^ < F.^.|V;:i-)(f,r.)^(-)(f,r^)T.[e^^ri/-''°'-g.^'}l]|F.g. > 

z<F.5.|T.[e^^--'^°"M7l]|F.5. > ^ ^ 

where [...] stands for the parametric time order. 

The computation of the single hole Green's function will be done in two steps: 

a) We compute first the parametric Green's function D^^^[t, K; g{T),T, A] where r is the 
correlation time interval and the coupling constant of the theory depends parametrically on 
the finite size T. For r > (holes excitations) the Green's function D^~^[t, K;g{T),T, A] 
is computed with the help of the effective action >5'f//~^" ''^[t j,, t/] defined on the finite time 



interval T. Using the method of finite size scaling [3J, l35|] we will compute the Green's 
function D^^^ [r, K; g{T),T, A] and the Fourier transform D^^^ [u, K; g(T),T, A] with respect 
the parametric time r at a fixed temporal size T and a fixed coupling constant g{T). 

b) The physical Green's function G^~^[T, K;g{T),T, A] is related to the parametric 
Green's function ^("^[r, K; g{f ),T, A]: 

G(-)[T,ie;^(f),T,A] = /o°° dT6iT-T)D(-^[T,K;gif),T,A] 

Once the Green's function D^^'^[uj, K; g{T),T, A] has been obtained, the physical Green's 
function G^^^[uj,K] is evaluated using the Fourier transform properties: 

2w 



G^-\u,K] = f dT{ / + / ) ^e^(— ')^D(-)[a;',i?;^(f ),T, A] (62) 
"'-f "'--^ •'jfi 
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Due to the finite size effect, the frequency integration oj' is restricted to £" > > 
where E is the band width and ^ is the finite size frequency cut-off. The bandwidth 
E = vpA is given hj vp = — x) where Kp is the Fermi momentum at half filhngs 

and X represents the hole doping. When the hole doping x increases, the bandwidth decreases 
E{x 1) 0. 

One of the interesting consequences of our formulation is that the single particle (hole) 
Green's function is a function of the effective time interval T. The coupling constant of 
the theory depends on the time interval between the creation and the destruction of the 
hole. Therefore we do not have one single action for all the time intervals. For example, the 
single particle Green's function for an infinite time interval ti — t f —>■ oo is described by the 
non interacting free action. Using the RG theory we compute the time dependent Green's 
function for a fixed time interval — tf = T. The frequency dependent Green's function is 
rather non trivial since we have to perform a time integration over all the time and over all 
the possible coupling constants! Performing the Fourier transform by integrating over all 
the time dependent coupling constants (which are a function of the time intervals T), we 
will show that the self energy is dominated at low frequencies by a relaxation part which 
is linear in frequency Tij^iuj) oc ujF{:^^). The function E{:^^) represents the crossover 
from 1 when u ^ to u for increasing frequencies. The crossover region is determined 
by the bandwidth function vpA. With increasing doping the bandwidth decreases and the 
crossover region shrinks to 0. As a result, S/m(co') is modified to the Fermi liquid behavior 
S/.m(u;) oc u"^. 



IX. THE RENORMALIZATION GROUP FOR THE ACTION 

S^ff-'^°^^[ti,ff;g{f),T, A] IN TWO DIMENSIONS 



The explicit form of the effective action in two space dimensions given in eq.(57) 
S'^-ff~'^°'-^[ti,tf;g{T),T, A] is restricted by the temporal size T. Therefore the method of 
finite size scaling will be used. 

The scaling dimensions for the coupling constant in eq.(57) are obtained from a mo- 



mentum cut-off which is normal to the Fermi Surface 



44j . The scaling dimensions in the 
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vicinity of the Fermi Surface is dp.s. = 1 (the scahng dimensions is only modified around the 
corners (0, ±7r) , (±7r, 0) close to half fillings). As a result, the coupling constant scales like 
g{T)/vp = g(T)E'^~'^''-^- = g{T)E. The fixed points for any theory are achieved by taking 
the limit 6 = e' ^ oo where 6 = describes the reduction of the bandwidth cut-off For the 
present problem the limit 6 = — oo can not be taken since we have to stop the scaling at 
b = bx = vpAT = t. 

The critical behavior g{T)E'^~'^^-^- is investigated using similar methods as employed for 
the I sing — gif^ model in d=3 dimensions. For the Ising case the coupling constant g obeys 
g = gE^~'^ = gE. (For the Ising — gip^ one performs the calculations at a fictitious 
dimension d = A — t] such that at the value rj = the coupling constant is marginal. The 
R.G. equations take the form ^ = (A — d)g... = 7]g...; to recover the physics for d=3 we take 
the limit rj —>■ 1 at the end of the calculation.) 

Following the analogy with the Ising — g(p'^ R.G. we introduce fictitious dimensions of the 
Fermi Surface dp.s. = 2— ?7, such that when rj ^ 1 one reproduces the one dimensional scaling 
of the Fermi Surface. The integration variable s] is replaced by ^^-^--n ^"i-r] ^ such 

that at the limit 77 ^ 1 we obtain J2-r,[e, s] — * J[e, s\. As a result, the scaling dimension of 
the coupling constant becomes marginal for 77 = 0. For < r/ < 1 we find: 



g{f) = ^(f )E2-(2-.) = ^(f (63) 
We will use the differential R.G. method where the integration in the energy shell 



E — dE < e < E is performed using the differential variable dl = ^ |32|, l33|, |43| . We find 
that the coupling constant constant g{T) for the zero angular momentum channel 43 1 
obeys the following R.G. equation. 



Comparing eq.(64) with the R.G. equation for singlet superconductivity [435, observe 
that due to the additional time integration in eq.(57), eq.(64) has a linear term g with the 
scaling dimension rj = 1 and that the term g"^ is rescaled by the temporal finite size parameter 
PiTE = t). The dimensionless parameter obeys 1 < (3{TE) < \T\E = \T\vfA = t. At the 
limit T — s> 00 one finds that the R.G. equation has an infrared stable fixed point g = 



An 



Vt3{TE) 



0, which describes the Marginal Fermi liquid. We use the complex representation 
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for the coupling constant g — —i{u + iA) with u — 3( = ^ and A = uSin{S), with 

the initial condition 5 — > 0. We find that the equations have an infrared stable fixed point 
given by {u* — > 0, A* = ^^^^^^^^ 0). Due to the finite temporal size T the R.G. equation 
(64) is only valid for < / < Log{\T\vFA) = Log{t) = k. 

In order to construct the full R.G. flow, we have to compute the differential self energy 
from which we will extract the wave function renormalization. We will work with the 
two dimensional representations of the coupling constant g — —iu + A, with the initial 
conditions A(Z = 0) — > 0. The self energy of our action is a function of frequency and 
couphng constants S[ci;, e; u, A, E\. We will use the infrared stable fixed point to compute 
the self energy S[c(j, e; -u, A, E] and to tune the chemical potential 5^f- We find that 5^p{s) 
is given by the same self energy at zero frequency for all the points s on the Fermi surface, 
5^f{s) — E[0, 0; M*, A*]. Expanding the self energy in powers of uj allows to compute the 
wave function renormalization Z^. We find: '^^"9lZj{u(i)] _ 

As a result, the previous R.G. equations for u and A are modified: 



du , u. A/ N uA 
dA , u. . ,A^ — 

- = i,--)A-m^^ (65) 

This set of equations have the initial conditions u{l = 0) = 3(^t^)^ = and A(/ = 0) = 
u{l = 0)Sin{6) 0. Due to the finite time interval T we have to restrict the scaling to the 
domain 1 < ^{TE = t) < TvpA = t. 

As a result, the coupling constant will reach the end point values u{lt) — u{Log{t)) and 
A{lt) — A{Log{t)). The finite size scahng results are given by the numerical solution of the 
R.G equations. 

This set of equations have an infrared stable fixed point given by {u* = 0, A* = 0). For 
a finite time interval T the values of the coupling constant deviate from the fixed point. For 
t — >■ CX3 we obtain that the coupling constants reach the values u{lt) = ^ and A{lt) ~ 
where ui and di are universal constants. The result of the R.G. fiow are given in figure 
1. In figure 1 we show that the dissipative coupling constant A{lt) fits the analytic form 
d{t) = 0.7588/t with the universal constant di = 0.7588. 
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d(t)=0.759/t 
0.0008- 




2000 4000 6000 8000 1000 



FIG. 1: (Color online)The Dissipative Coupling Constant A{lt) : We show two graphs one is the 
R.G. result and the second is the fit to the analytic form d{t) = 0.7588/t = di/t 

X. COMPUTATION OF THE PARAMETRIC GREEN'S FUNCTION D FOR 
TIMES T<T 



Using the R.G. results from the previous section we will compute the hole type parametric 
Greens function for parallel spins D^~^[uj, e]g = —i{u + iA)] = D^~'^[uj, e]u, A], where e is 
the energy excitation perpendicular to the Fermi Surface e = vp ■ {K — Kf{s)). 

This Greens function [cj, e; m, A] is computed using the unperturbed Green's function 

{D^-\u, e; u, A, E])-' = {D^^ [u, e])-' + S[^, e; u, A, E] 

= LU -e-5fiF + J:[0,0; u, A, E] ^' (66) 

duj 

where Ti[uj,e;u, A, E] is the self energy. Using the fact that the R.G. equation has an 
infrared fixed point we can tune the shift in the chemical potential 6fiF{s) such that 6fiF{s) = 
S[0, 0; u{lt) ^u*^f, A{k) ^A*^^,E]. 

We obtain the wave function renormalization Z^(u{lt)) = e~ -fo ^'^^ . In the present case 
we stop scaling at = vpAT = t and we find: 

Z^uil = log{t) = h)] = e" Jo )'^' ^ e"^ (67) 

As a result we obtain the finite size Green's function for t » 1 with the universal 
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parameters ui and di: D'^ ^uj.e-.uilt), /S.(lt),E\ = — ^ "^'a^.i ~ — - — ^—dr- The action in 
eq.(57) is restricted for large time intervals T > for which we can replace e ~ 1. 

The parametric Green's function which is restricted to the frequency interval < uj < 
vpA is given by: 

D[u;, e; u{k), A{k), E] = + (68) 

ij — e + to cj — e + lOJ^ 



XI. COMPUTATION OF THE PHYSICAL GREEN'S FUNCTION G 



The Physical Green's function G will be computed from the Fourier representation of the 
parametric Green's function D given by equation (68). We substitute in equation (62) the 
explicit form of the parametric Green's function D[w, e] = D^eiw, e] +iDim[w, e] as given in 
equation (68) . 

Using the Fourier transform of equation 62 we find that the Green's function G is given 
in terms of two new functions Fji{t;u!,e) and Fj^{t]uj,e) which are a linear combination 
of the real and imaginary part of the parametric Green's function D'^'~\w^ e] = D^^^\w^ e] + 

G^Rj[w,e]= J^FR{t-u,e)dt 

vpA 

G^j-J[w,e] = J^FjUt;u,e)dt 

vpA 

where 



FR{t;u;,e) ^ —[D\J[zu,e]Cos{{u;-u;)t)-D\J[-zu,e\Cos{{u; + zu)t)] 

rv. 



-[D'fJiw, e]Sin{{u; - w)t) - D^-J[-w, e]Sin{{u; + w)t)] (69) 



rVpA , , , , 

Fim{t;uj,e) = —[D)iJ[w,e]Sin{{u-w)t)-D)iJ[-w,e]Sin{{uj + w)t)] 
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- — [Dt^Ke]Co.((^-a;)t) -Dt^[-tu,e]C'o.((a^ + t^)t)] (70) 

The Physical Green's function G[uj, e = vp ■ {K — Kf{s)); E = vpA] is given in terms of 
the self energy T,{uj, e) = S/je(ci;, e) + iT,jrn{^, e)- 

G[uj,e,E - — H — — — ^ (71 j 

In the present case equations (69) and (70) gives us the Physical Green's in terms of 
the parametric Green's function D. We represent the self energies of G in terms of the 
parametric Green's function. We find that: 

ERe{uj,e; z) = GRe[w,e; z]/ {{GRe[w, e; z]Y + {Gjm[w, e; z]f) - w 
and 

S/^(u;,e; z) = -Gim[w, e; z]/{{GRe[w, e; z]y + (G im[w , e; z]y) - z where z ^0 
where GreIw, e] and Gimlw, e] are given by the equations (69) — (70) . 
The results for the self energy are given in figures 2 and 3 at a fixed energy e = 0. 
(The Green's function will be given as a function of dimensionless frequency and energy 
gj;^^^ uj and ^^^^^ e.) We observe that with increasing frequency the calculation 
of the self energy at energy e = becomes less accurate. For larger frequencies we have to 
compute YjReiy^^e) and S/m(^^,e) at finite energies e. 

The self energy S/m(u;) oc loF{^-^^ is linear in frequency for ^ 0. The function F{^-^^ 
represents the crossover from 1 when ^ to for increasing frequencies. When doping 
increases, the bandwidth vpA decreases and we observe that the linear region in frequency 
shrinks to 0. As a result HiyniiS) is modified to the Fermi liquid behavior S/m(ti;) oc lJ^ . 

Figure 3 shows that the imaginary part of the self energy for low holes densities is linear 
in frequency. As a result, the single hole excitation (at low frequencies) has a width which 
is linear in frequency, and the scattering rate obeys ^ oc u;F(^j^) in agreement with the 



infrared data 25 1. 
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FIG. 2: (Color online)The real part of the self energy self energy of the physical Green's function 
G{uj, e) computed from the finite size Green's function e; u{lt), A(/i), E) 
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FIG. 3: (Color online)The Imaginary part of the self energy of the physical Green's function G(a;, e) 
computed from the finite size Green's function D{uj, e; u(/t), A(/t), E) 

XII. APPLICATION OF THE THEORY TO SUPERCONDUCTIVITY 



In this section we will attempt to connect the theory presented with the physics of the 
high Tc material. In particular we have considered a model at zero temperature away from 
half fillings where the magnetic order has vanished. Under this condition we have shown 
that the marginal Fermi liquid with an imaginary self energy which is linear in frequency 
is obtained. This results are in agreement with the experiments which show that at T=0 a 
window exists between the magnetic ordered state and the appearance of superconductivity. 
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Therefore, when the exchange interaction is included a d-wave superconducting [26| phase 
within the exclusion of double occupancy will appear. The d — wave order parameter 
I]a['^T(^)'^i(^ + a) — ip^{x)ipi^{x + a)]|\E'° >^ {contrarily to the s — wave) is compatible 
with constraints ■?/'|(x)^/'|(a;)|\l'° >= 0. Further doping of the superconductor at T=0 will 
give rise to a transition from a free vortex monopole phase to a spin wave phase jsj. 

For the remaining part of this section we will show that the effect of exclusion of double oc- 
cupancy gives rise in the superconducting phase to an asymmetry in the tunneling density of 
states. In order to demonstrate the asymmetry effect of the projection of double occupancy 
we consider a qualitative calculation for superconductors . Strictly speaking an accurate 
comparison with the experiment must use the full d wave structure. In order to demonstrate 
the effect of asymmetry induced by the projection, it is important to show that the asymme- 
try can be obtained also for an uniform state. We consider the standard BCS hamiltonian 
and use the projection introduced in the previous sections. In the absence of the projection 
the effect of the BCS gap Abcs gives rise (after integration over the single particle energy) 
to the following tunneling density of states, N^^~'^"'^'^\uj) = f / dtlm.G^^^ {t,uj)\ 



^(5-«,a.e)^^^ = f' de[{l + ^ ^ )(- 



+ (1+ , )( , )] (72) 



where z ^ 0. 



Next we repeat the calculation when we project out double occupancy! For this purpose 
we use the self energy given in figure 3, S/m(u;). 

The tunneling density of states for the projected case of double occupancy is given by 



Nt''-'-\uj) = f de[{l + , ^ )(- 



y2 



+ (1+ , ' )( ^ )] (73) 



where z ^ 0. 
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FIG. 4: (Color online)a)-The tunneling density of states for S-wave superconductivity as a func- 
tion of the voltage V without projection Nj? u)ave)^eV^ b)-The tunneling density of states for 
S-wave superconductivity as a function of the voltage V with projection (the asymmetric graph) 



{S-w.ex.)feV\ 
V h > 



In order to emphasize the asymmetric effect of the self energy we consider typical values 
of temperatures and gap. For the gap we take the value Abcs = 0.38 x 10^'^ eV and restrict 
the temperature to, Tb < ^bcs- 

In figure 4 we show the two graphs of the tunneling density of states, N^~'^"'^^\uj = ^) 
is the tunneling density of states in the absence of projection and A^^'^ ^'^^'\uj = ^) is 
the tunneling density of states for the projected case. The tunneling density of states 
as a function of the tunneling voltage V shows a clear asymmetry between the projected 



(S—w.ex.) 



and the non - projected function A^^ 



(S—wave) 



In figure 5 we show on the same graph: the experimental data for the tunneling of 
the density states observed in {27] and our projected tunneling density of states. Figure 5 
demonstrates that the asymmetry in the tunneling density of states can be explained by the 
projected Green's function. 
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FIG. 5: (Color oiiline)The tunneling density of states given by the dots extracted from the exper- 
iment in ref. I27I compared to the analytic formula for the projected tunneling density of states 

XIII. CONCLUSION 



We have proposed a new method to study correlated electrons where the traditional 
method of slave particles is avoided. We have demonstrated that by choosing proper variables 
additional constraints can be included. We obtain a system of first class constraints which 
generate gauge transformations. 

The authors [l^ have computed the wave function using only one constraint. We 
identify two additional constraints which form an non-Abelian group and neglect according 
to the R.G. analysis the secondary constraints. As a result the simple delta function con- 
straint is replaced by a non-linear integration measure which gives rise to time dependent 
action. 

The effective action is analyzed with the help of the R.G. method for finite size systems 
(in the time domain). 

The R.G. analysis allows to compute the single particle self energy which is used for 
qualitative comparison with the experiments. 
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APPENDIXiTHE SECONDARY FIRST CLASS CONSTRAINTS 



In this Appendix we will show that the effective interactions induced by the secondary 
first class constraints are irrelevant operators for describing the long low energy Physics and 
therefore can be neglected. 

In order to show this we perform the foUowings steps: 

Al) Compute the commutator of the kinetic energy with the primary first class con- 
straints. (In the absence of the exchange interaction, the hamiltonian is given by Hq.) 

The secondary constraints q{x), q^{x) and q^ix) are obtained by commuting the primary 
constraints with the hamiltonian Hq and subtracting the primary constraints: 

q{x) = ][Q{x), Ho] - Q{x) = Ea['^a=l{x)ilJa=]{x + 3)+ ^a=l{{x + a))^^=|(f)] - Q{x) 

q+{x) = |[Q+(f),//o]-g+(f) =Ea[V'a=i(^)V'.=T(^+«)+V'a=i((^+«))V'a=T(^)]+-Q+(^) 
and 

qs{x) = ^[Qs{x),Ho] - [Qsix) - 1)] = Ea['0a=i(^)V'a=i(^+ «) + '0a=T((^+ «))V'a=T(^)] - 

[QsiS) - 1] 

A2) In section VII we have parametrized Fermi-Surface in terms of the polar angle s, 
normal N{s) to the Fermi Surface and chiral fermions Rfj{x,s), ,Li^[x,s) .The hamiltonian 
Hq is given by: 

^0 = Id^xJ^ ^ Ea=tM{x; s){-tVFh)N{s)-d^K{x; s)+LUx; s){tVFK)N{s)-d^L,{x- s)] 
A3) Approximating the difference between the commutators and the primary constraints 
by a first order spatial derivative around the Fermi-Surface we obtain: 

q{x) ^ ^[R^{x, s)N{s) ■ dgL^ix, s) - Li{x, s)N{s) ■ dgR^{x, s)] 

Q^i^) ~ [Jo f [^t(^, ^)N{s) ■ dgL^{x, s) - Li(f , s)N{s) ■ d,R^{x, s)]]+ 
and 

qsix) ^ ^[Rtix, s)N{s) ■ d^R^ix, s) + s)N{s) ■ d^L^{x, s)]. 

A4) The presence of the secondary constraints the modifies the Lagrangian L ^ L + SL 
where SL is given by: 
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5L = E'=i 7r.(f, t)dty{x, t) - E'=i(7r.(f , t)r{x, t) + y{x, t)q(-\x, t)) 
where X'^{x,t) are the Lagrange multiphers which enforces the secondary constraints, 
TTr^x, t) are the canonical momentum conjugate to the new Lagrange V'(a:, t) multiphers and 
qlr\x,t)) arc obtained using the linear transformation given in equations (35) and (36). 
This gives rise to the evolution operator: 

In order to compute the effective interaction induced by the secondary constraints we 
need the integration measure for the secondary Lagrange multipliers X''{x.t). We approxi- 
mate the measure by a regular integration and obtain a set of delta functions which enforces 
the constraints q^~\x,t)). The delta functions constraints effectively replaced by expo- 
nentials of Gaussian terms with the coupling constant Sgj { which at short distances goes 
to infinity and therefore is equivalent to a delta function). The Gaussian action is given 
by; 5giY.f=iI d'^^ I dtij dt2[qi~\x,ti))ql~\x,t2))]. As a result, we obtain the correction 
SS^ll~'^"'''^[ti,tf] to the original action given in eq.(57): 

^Siil~''"^%,tf] = I — I — U dti [J dt2 n i^r\3Jh,Si]J[e2,S2] 

Jo IT Jo IT Jti Jti [ZTT)-^ 

5gi5{-ei + 62 + 63 - 64) 

"^(ti, -ei; si){e2)Lf'\h, 62] Si)4"^(t2, 63; S2)i-ei)R[-\t2, -64; S2) 
+L\^~\ti, ei; Si){-e2)Rf~\h, -62; Si)L\-\t2, 63; S2){-e^)R{-) ^{t2, -64] S2)] 
+[Rf'\h, -ei; Si){-ei)Lf-\h, 62; Si)R\-\t2, -63; S2) (64)4-^(^2, €4; S2) + 
Lf~\tu ei; si)(ei)i?f '^(ti, -es; Si)L\-\t2, e^; S2) (-64)4^^(^2, -e^, S2)]].... 

The spatial derivatives in the last equation are replaced by the energy excitations normal 
to the Fermi Surface e. As a result, the engineering dimensions of the coupling constants Sg 
is given by: Sgi — SgiVpE'^. 

The presence of two spatial derivatives and the two time integrations generate the engi- 
neering dimensions E~^. Therefore we will ignore SSlH"^"^^ for describing the Physics at 
low energies. (The Physics at short distances is sensitive to operators which have negative 
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scaling dimensions and therefore can not be ignored.) 
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